Tracer diffusivity in a time or space dependent temperature field 

Ramin Golestanian 1 ' 2,3 and Armand Ajdari 2 
1 Institute for Advanced Studies in Basic Sciences, Zanjan, 45195-159, Iran 
2 Laboratoire de Physico-Chimie Theorique, UMR CNRS 7083, E.S.P.C.I., 75231 Paris Cedex 05, France 
3 Institute for Studies in Theoretical Physics and Mathematics, P.O. Box 19395-5531, Tehran, Iran 

(to appear in Europhys. Lett.) 

The conventional assumption that the self-diffusion coefficient of a small tracer can be obtained 
by a local and instantaneous application of Einstein's relation in a temperature field with spatial 
and temporal heterogeneity is revisited. It is shown that hydrodynamic fluctuations contribute to 
the self-diffusion tensor in a universal way, i.e. independent of the size and shape of the tracer. 
The hydrodynamic contribution is anisotropic — it reflects knowledge of the global anisotropy in 
the temperature profile, leading to anisotropic self-diffusion tensor for a spherical tracer. It is also 
retarded — it creates memory effects during the diffusion process due to hydrodynamic interactions. 
PACS: 05.40.-a, 66.10.Cb, 05.60.Cd 



Introduction - Motion of particles in thermal gradi- 
ents is a relatively old subject, of both fundamental and 
applied interest. The net directed motion in a steady 
gradient (usually referred to as the Soret effect or ther- 
mophoresis) has been observed and quantified for vari- 
ous systems (polymer molecules, solid particles or bub- 
bles in a fluid) B but often through indirect measure- 
ments. However, the microscopic mechanisms responsi- 
ble for this motion are far from clear in many cases (see 
e.g. Ref. Q). The resulting collective motion is usually 
rationalized through Onsager-like relations, with kinetic 
coefficients that are often empirically deduced from indi- 
rect measurements. The corresponding experiments are 
often difficult to set up due, in particular, to the buoy- 
ancy driven convective effects induced by thermal gradi- 
ents. 

In principle, the development of microfluidics (see H^|] 
and references therein) should allow to create fluidic envi- 
ronment less prone to such convective instabilities, where 
temperature gradients can be induced and controlled, 
and that should allow direct measurement of particle mo- 
tion. Various fluorescence techniques have indeed been 
developed that allow a rather precise detection of temper- 
ature (|] , and a laser in the IR range can offer a versatile 
way of inducing these gradients. 

Although the directed motion of particles (Soret effect) 
is the most obvious target of such a study, its amplitude 
often scales with a high power of the size a of the particle 
(typically a 3 ), so that for small tracers an equally inter- 
esting question is that of the diffusivity of the particle. 
In this paper, we take a first step in this direction on 
the theory side, and use the simplest model possible (an 
incompressible liquid in the Stokes limit with thermal 
fluctuations of the velocity field, where the temperature 
profile is externally imposed and not modified by ve- 
locity fluctuations, and where viscosity and density are 
assumed to be temperature independent), to study the 
effect of spatial and temporal variations of the tempera- 



ture on the diffusivity of a tracer 



Main results - More precisely, we re-examine the as- 
sumption of local and instantaneous dependence of dif- 
fusion coefficient on temperature, by taking into account 
the effect of hydrodynamic fluctuations at low Reynolds 
number. We have studied the diffusion of an arbitrarily 
shaped particle (tracer) with mobility tensor flij in a fluid 
with mass density p and viscosity r/, which undergoes hy- 
drodynamic fluctuations that are driven by a given time 
and space dependent temperature profile T(r,i). (For 
example, fi%j = 8ij / {Qitrja) for a spherical particle of ra- 
dius a, while it can be a more complicated anisotropic 
tensor for an arbitrarily shaped particle.) We find that 
the Einstein formula for the self-diffusion tensor takes on 
the form 



Aj(r, t) = ti tJ k B T(r, t) + SD^r, t), 



(1) 



with a (usually small) correction term that is univer- 
sal, i.e. independent of the shape of the tracer particle 
[gj. This term is caused by an interplay between the 
long-ranged hydrodynamic interactions and the spatial 
or temporal heterogeneity of the temperature field. For 
a steady space dependent temperature profile T(r), we 
find 
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which implies that a global anisotropy in the temperature 
profile will be reflected in the self-diffusion of a tracer. 
For the case of a temporally varying temperature T(t), 
on the other hand, we obtain 
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which manifests a memory effect corresponding to retar- 
dation in hydrodynamic flows. 

Model - We consider an incompressible fluid at low 
Reynolds number that is subject to a thermal random 
stirring force (density) f(r, t). The stochastic dynamics 
of the velocity field v(r, t), which is subject to the incom- 
pressibility constraint V • v = 0, is then governed by the 
linearized Navier-Stokes equation 



pdtv = —Vp ■ 



rjV 2 v 
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in which p is the pressure field, and the density p and vis- 
cosity rj are taken to be temperature independent. The 
noise in the above equation is assumed to be Gaussian, 
with (fi(r,t)) = and 

{f^v^f^v'.t')) = 2r,fcBnr,i)M-V 2 )5 3 (r - - 0. 

(5) 

where the noise correlator is chosen so as to ensure local 
and instantaneous thermal equilibrium for a sufficiently 
slowly varying temperature field Q . Note that the struc- 
ture of the noise term should be such that it only induces 
vorticity-free velocity fluctuations, so that the system is 
not driven far from equilibrium by the noise |10|. The 
velocity field correlation function can then be calculated 
in a straightforward way. We obtain 

(«»(!■, *)« i (r / ,i / )) =2r/ / dh d'riAfelXri,*!) 



x VG fei (r - n,i - h) ■ VG kj (r' -r 1? t'-tx), (6) 

where GV,- is the Green's function for Eq.(|]) above, i.e. 
in Fourier space Gy(q, u>) — (<% — qiqj)/(—ipuj + rjq 2 ). 

Using the result for the velocity correlation function, 
we can extract the diffusion coefficient of a tracer by as- 
suming that the velocity fluctuations of an immersed par- 
ticle would follow that of the fluid at the particle's loca- 
tion by way of no-slip boundary condition. This assump- 
tion is reasonable as long as we are probing the long time 
behavior of the motion. We thus define the self-diffusion 
tensor via 



dt' (« i (r,t)« J (r,* / )), 



(7) 



where it should be clear that this Eulerian velocity auto- 
correlation function decays rapidly. 

Examining the expression in Eq.(^) then shows that 
the diffusion tensor, as defined above, is decomposed as 
in Eq.(0), in which p^ = J ^rpp(% - hk^/^k 2 ) in- 
volves a diverging integral that is sensitive to the details 
at short length scales, whereas 5Dij is a long-ranged 
convolution of a hydrodynamic response function and 
the temperature field. Note that SDij = for a constant 
temperature profile. One can easily see that in fact pij 



corresponds to the mobility tensor of the particle, which 
we would obtain if we treated the proper hydrodynamic 
problem of flow past a body, as the form of the integral 
expression suggests if we use a cutoff of the order of the 
size of the particle In short, we can say that the 

diffusivity tensor Dij(r,t) obtained from Eqs.(||) and (Q) 
consists of two contributions: (1) a term that is sensitive 
to short length-scale features and thus to the geometry 
and size of the tracer, which depends only locally on tem- 
perature, and (2) a term that is sensitive to the non-local 
features in the temperature profile, and is thus indepen- 
dent of the size and shape of the tracer as long as these 
features have characteristic length-scales that are much 
larger than the size of the tracer. While the closed form 
expression for SDij(r, t) can be obtained for an arbitrary 
temperature profile, we choose to report here separately 
the results for space dependent temperature as in Eq. (||) , 
and for time dependent temperature as in Eq.(|J) above. 

Steady non-uniform temperature - Let us now try to 
elaborate on the meaning of the above results. For a 
non-uniform temperature, Eq.(||) implies that a spheri- 
cal particle should undergo anisotropic diffusive motion, 
with the anisotropy being inherited from the structure of 
the temperature profile For example, a harmonically 
modulated temperature profile T(r) = Tq + 5Tcos(k • r), 
yields SD^r) = -fc B <5Tfc/(64?7) (3% + fc i fc i )cos(k • r) 
that predicts an anisotropy of SDn/SD± = 4/3 for the 
ratio between the corrections to the diffusion constant in 
the parallel and perpendicular directions respectively. It 
also shows that for a small tracer of size a the correction 
is typically weaker than the usual term in (1) by a fac- 
tor of order ~ (ka)(\8T\/To), with both terms typically 
smaller than 1. 

To make connection with experiments, we need to con- 
sider realistic temperature profiles. Non-uniform temper- 
ature can be achieved by locally heating the fluid, and 
the profile can be obtained through the heat equation 
— KV 2 T(r) = Q(r), where k is the thermal conductivity 
and Q(r) is the volume density of heat generated by ex- 
ternal sources per unit time. The heating can originate 
from Ohmic loss if, for example, a wire carrying elec- 
tric current is placed in the fluid, or from optical heat- 
ing with an IR laser light focused into a spot within the 
fluid. Note, however, that in the latter case we need to 
make sure that the dielectric constant of the tracer par- 
ticle matches that of the fluid at the wavelength of the 
beam, or that the particle is sufficiently far from the spot, 
so that optical trapping of the dielectric particle does 
not interfere with the effect that we are studying here. 
We choose a source with the anisotropic Gaussian pro- 
file Q(r) = P /(£ 2 £ z )exp{-Tr{x 2 + y 2 )/£ 2 ~irz 2 /£ 2 z }, in 
which Pq is the total dissipation power, and find a va- 
riety of different asymptotic regimes: (i) For £,£ z <C r, 
where the tracer is sufficiently far from a point-like heat- 
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ing source, we find T(r) = Tq + Po/(4:irKr), and 
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to the leading order. Note that SD has a different (faster) 
power law decay with r as compared to that of T. More- 
over, we find an anisotropy SD rr /SD±± = 1/2. (ii) For 
I <C r <C £ z , where we have a linear source, we find 
T(rj_) = To — (Po/2wk£ z ) ln|rj_|, and corrections to the 
diffusion tensor as 



SD zz (r ± ) = - 
5D±{t x ) = 
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and 5D z ±(r±) = to the leading order. Since in this 
case the temperature variations in space is very weak 
compared to the variations in SD, it may be a good can- 
didate for experimental observation of the space depen- 
dency of SD. We also find an anisotropy SD zz /SDge = 
SDpp/SDeg — 3/4 in terms of the components in cylin- 



drical coordinates, (iii) For t z 
a planar source, we find T{z) 
the corrections as 
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and 5D z ±(z) = to the leading order. This case also 
corresponds to a temperature profile with a constant gra- 
dient, and is expected to have the strongest effects in 
terms of magnitude as compared to other regimes. The 
anisotropy in this case is found as SD ZZ /SD±± = 4/3, 
similar to the harmonically modulated case, (iv) Finally, 
for I, £ z 3> r, where we have an extended anisotropic 
source with an aspect ratio a = £ z /£, we find T(r) = 

Q+\/a 2 — 



T + (P /2ttk£) In 
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to the leading order, 
3a 2 )tanh _1 (-s/i - a 2 )~VT 



where A z = [(4 — 

^] /(l-a 2 ) 3 / 2 , and A ± - 

[(6 - 7a 2 )tanh- 1 (VT _ + Vl^?)/[2{1 - a 2 ) 3 / 2 }. 
This is a very interesting case in the sense that both 
temperature and diffusion tensor have constant profiles 
to the leading order, and yet the diffusing tracer feels 
the large-scale anisotropy in the temperature profile, as 
manifest in Eq.(|ll"|). The anisotropy in this case is found 
as 
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which depends on the aspect ratio, as plotted in Fig. 1. 
Note that Eq. ( |I2| ) shows that the anisotropy asymptotes 
to 4/3 in the limit £ z <C £, similar to the planar and the 
harmonically modulated cases. Moreover, it reveals that 
oblate profiles lead to stronger anisotropy than the pro- 
late cases. This trend is in agreement with diminished 
anisotropy in the case of linear sources {SD zz /SD pp = 1) 
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FIG. 1. A plot of the anisotropy in the correction to 
the diffusion coefficients Eq.(12) as a function the aspect 
ratio of the heating profile. The anisotropy is more pro- 
nounced for oblate profiles as compared with the prolate 
cases. 

Temporally varying uniform temperature - For time de- 
pendent temperatures, hydrodynamic retardation leads 
to non- locality in time, as evident in Eq.(||) above. This 
memory effect can be well illustrated in the case of an 
oscillating temperature field of the form T(r) — Tq + 
ST cos wot, where it leads to an isotropic correction to 
the diffusion tensor as SDij(t) = 8ijSD(t) with 
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The above result shows a systematic phase lag of 37r/4 for 
the hydrodynamic correction to the diffusion constant, 
with respect to the temperature. The characteristic 
length in the above equation is set by the kinematic 
viscosity rj/p (that has the dimension of a diffusion con- 
stant) and the frequency ujq as ot = yV/Mo, and is of 
the order of 1.3 mm for water at luq/(2tt) = 0.1 Hz. 

Discussion - The above results are interesting from a 
fundamental point of view because they reveal that hy- 
drodynamic interactions reflect on a diffusing particle the 
global knowledge of the temperature heterogeneity and 
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anisotropy. They may appear though to be difficult to 
verify experimentally, because in typical cases the correc- 
tion will be considerably smaller than the conventional 
value given by Einstein's relation. However, we would 
like to point out that due to the anisotropy in the dif- 
fusion tensor, one can choose to probe a quantity such 
as D xy , by direct measurements of (Ax(t)Ay(t))/2t, in a 
non-symmetrical coordinate system. Since this quantity 
is always zero for a spherical tracer in the local picture, it 
would provide a feasible experimental test of this effect. 

The structure of the stochastic forcing in Eqs.(^) and 
(||) has a crucial role in determining the form of the above 
results. While we are studying systems that are only 
slightly away from equilibrium (by choosing a noise cor- 
relator that is a smooth perturbation of the one that en- 
sures thermal equilibrium), other noise structures could 
drastically change the results so that for example the de- 
composition in Eq.(jl|) could be blurred and the scaling 
with length could be altered. 

Let us recall that our description of the fluid dynamics 
is clearly a great simplification, which could be improved 
by taking into account the temperature dependence of 
the density (and thus release of the incompressibility as- 
sumption) and the viscosity, and also the convection of 
thermal fluctuations, along the lines of the approach used 
in Ref. |l3| to study fluctuations of a fluid away from 
equilibrium. With these ingredients, the long-range char- 
acter of the effects described here should then suffer from 
screening, but the corresponding screening length will re- 
main large for weakly compressible fluids. 

The hydrodynamic contribution to the effective diffu- 
sivity is independent of the size and shape of the tracer 
as long as the shortest length-scale over which the tem- 
perature profile varies significantly is considerably larger 
than the size of the tracer particle. In this regard, a 
proper account of the tracer geometry for larger parti- 
cles in strong temperature gradients may become impor- 
tant, when the interplay with Soret effect is expected to 
appear more clearly. We nevertheless expect that some 
qualitative features obtained in the present work should 
remain, i.e. that the fluctuating character of the particle 
motion should be sensitive to the geometry and history 
of the temperature field. 

Eventually, as pointed out in the introduction, the 
most important steps to be taken are likely to be on the 
experimental side with either direct observation of single 
tracers, of sets of particles in a well-controlled fluidic and 
thermal environment, or direct probing of the anisotropy 
in the solvent self-diffusivity (in which case we expect 
minimal interference due to Soret effect) Q. 
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